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Abstract : The models of translationally invariant infinite nuclear matter in the rela- 
tivistic mean field models are very interesting and simple, since the nucleon can connect 
only to a constant vector and scalar meson field. Can one connect these to the compli- 
cated phase transitions of QCD ? For an affirmative answer to this question, one must 
consider models where the coupling constants to the scalar and vector fields must depend 
on density in a non-linear way, since as such the models are not explicitly chirally invari- 
ant. Once this is ensured, indeed one can derive a quark condensate indirectly from the 
energy density of nuclear matter which goes to zero at large density and temperature. 
The change to zero condensate indicates a smooth phase transition. 
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1. Introduction 

There is lot of interest in deriving QCD parameters from various relativistic hadronic 
model calculations [|TJ, 0, [[|. In the present paper : 

(1) we highlight the fact that different hadronic models, - which fit the binding energy 
S/p — M at the density of the infinite nuclear matter p Q , - also gives us the same quark 
condensate at p Q . But the expected behaviour, of (qq) — > at high p and T, is not built in 
the models. Indeed since the models are not explicitly chirally invariant, one can combine 
the behaviour (qq) — > with an effective nucleon mass which is not changing very much 
at high density. In the next section we will discuss this in more detail with respect to the 
recent queries |4[] regarding chiral symmetry restoration. 

(2) qualitative arguments are given to relate the behaviour of quark condensate (qq), 
to the incompressibility of nuclear matter in hadronic models (HM, in short). 

(3) we stress that to link QCD to hadronic models is important. The current at- 
tempts to derive QCD parameters from HM may provide important clues that will help 
in derivations of NN force from QCD. In QCD there is only one parameter Aqcd in terms 
of which one gets the strong coupling a s . This in turn must control the quark and gluon 
condensates, although to our knowledge can only cite heuristic arguments to show this, 
for example in ref. B : 

m dyn = Se^Agco = 300 MeV (1) 
where m dyn is the consituent quark mass and 
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-na s (qq) Q = -(m dyn f (2) 

with A Q cd = 130 MeV. 

In nuclear matter, there are more parameters, the binding energy -15.75, the equi- 
librium density and the incompressibility which are now leading to a specific value of 
(qq) p , replacing them, in a sense, in terms of one parameter. The hadronic models which 
we consider, namely the Walecka model |5J and its modifications |7[], show saturation of 
infinite nuclear matter through large cancellations between the attractive scalar and re- 
pulsive vector potentials. However this underlying NN force is as yet unknown from the 
QCD point of view, and in some way the QCD - construction of the NN force may be 
helped by the connection between (qq) and the incompressibility. 

The plan of the paper is as follows : in section 2 we discuss the quark condensate, 
how it is supposed to relate the scaling of the nucleon effective mass in a medium, the 
Hellmann-Feynman way to evaluate it from the quark condensate. In section 3 the nuclear 
models are discussed. Section 4 is devoted to results and discussions at finite density while 
section 5 discusses the nuclear equation of state. Section 6 contains the finite T results 
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for zero baryonic density but non-zero scalar density, which we believe is and will be 
probed in mid-rapidity reactions with very heavy ions. Section 7 contains the summary 
and conclusions. 



2. The quark condensate 

The quark condensate in QCD is : 

(qq) Q = -(230 ± 30 MeVf (3) 

The non-zero value is due to the breaking of approximate chiral symmetry by the vacuum, 
otherwise enjoyed by the QCD Hamiltonian, by virtue of the smallness of the quark mass 
m q : 

Hqcd = H + 2m q qq (4) 

the major part of the above being the chirally symmetric Hq. 

A necessary but not sufficient condition for chiral symmetry restoration in QCD is 
that the (qq) — > at high p and T |4]]. If the symmetry is still not restored the pseudo- 
Goldstone particle, the pion, will be massless and one can have pion condensation ||. 
Stretching the point further, can one presumably say that the effective nucleon will also 
be very light, if the pion becomes light ? Keeping these in mind we study the density 
and temperature dependence of the condensate and the corresponding sigma term. Let us 
begin with medium effects. Lattice studies cannot cope with finite density problems yet 
and QCD sumrule determination of the nucleon sigma term is not satisfactory- according 
to ref . . In this work hadronic models of nuclear matter are employed for our purpose. 

In the linear Walecka model ||, the condensate goes down linearly upto 1.5 times 
the nuclear matter density but then tends to increase, contrary to expectations based on 
ideas of chiral symmetry restoration. This is the same kind of result obtained by Li and 
Ko || very recently from Bonn potential in the relativistic Dirac-Brueckner approach. 
They also find the results unpalatable and conclude that effects like dependence of meson 
-nucleon coupling constants on the current quark mass and chiral invariance may become 
crucial in obtaining a reliable result for the density dependence of the quark condensate. 
Indeed we find the contrary result for variants of the Walecka model given by Zimanyi and 
Moszkowski (see 0), the models are labelled by ZM and ZM3 respectively where mesons 
interact non-linearly and couplings are p-dependent. In ZM3, the effective nucleon remains 
massive but (qq) goes to zero at high p. 

Before details of our calculations we discuss a few relevant points. In the QCD sum rule 



approach, studied by Ioffe JTU| and others [Tl[] one gets a simple approximate expression 
for the nucleon mass : 

M N = -(8n 2 /M 2 )(qq) (5) 
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and the formula is to be evaluated for M 2 ~ M^ 2 . This shows that to this order of 
approximation ( about 10 per cent or so) the nucleon mass is controlled by the quantity 
in eq.(|3]). From this rough analysis one would conclude that at finite temperature T or 
density p one should have a scaling 



[M N */M N ] 3 (6) 



where M^* is the effective nucleon mass in finite p or T. Indeed this was the idea put 
forward in [ F2| . Later it was suggested [I] that more careful analysis of the sum rule yields 



a linear scaling : 

(an) m 

[M N */M N ] (7) 



(qq) 

at least for finite density. In Figs. la and lb neither eq.(^|) nor eq.(^) is strictly valid 
beyond normal nuclear matter density po and the cubic scaling is roughly obeyed by the 
ZM model (Fig. la) while the linear one is preferred by the ZM3 model (Fig. lb). We 
now review how these figures can be obtained from standard nuclear matter calculations. 

A concrete way of getting the quark condensate was laid down by using the 
Hellmann-Feynman theorem : 

(^(rng)\^H Q c D \ip(rng)) = -^—(4){m q )\H Q c D \^{m q )) (8) 

which on using eq.(|4]) yields 

2m q {{qq) p - (qq) Q ) = m q —. (9) 

where the subscript p and indicate expectation value of the relevant operator for the 
nucleon in nuclear matter with uniform density p and for the vacuum respectively. The 
expression is the same for finite T. 8 is the appropriate energy density. The leading term 
of the expression d£/dm q of eq.(Q) is the experimentally known quantity <Ttv @ : 

2m q J d 3 x((N\qq\N) - (qq) Q ) = a N = m,-^. (10) 

where \N) is the free nucleon at rest. In nuclear matter of volume V with A = pV 
nucleons, translational invariance makes the quark condensate density constant and the 
integral in eq.flTOD gives : 

a A = 2m q V((qq) p -(qq) ). (11) 
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Comparing eq.(|TT|) to eq.(|9[) : 

p d£ 
A q dm. 



Oa = m q - — (12) 



"q 

and (using the relation of Gell-Mann, Oakes and Renner) : 



(qq) Amjtf P m ^tf 



with f w = 93 MeV and 



a A d(£/p) d(£/p) 

a °" = t = m "^r = aN ^M7 (14) 



is the effective a-commutator for a nucleon in the nuclear medium. If we neglect 5£, the 
nucleon kinetic and interaction energy density in 

£ = P M N + 5£ (15) 

we get a ef f = a N . 

Combining the eqns. (§, |ll| and [15]) we separate the static nucleon part and the effect 
of S£ on a e ff is : 

<?eff = <?n - -7— -r • (16) 



3. Relativistic nuclear matter models 



To consider in- medium quark condensate, one needs relativistic nuclear matter p[. 
Serot and Walecka advocated a relativistic field theoretic description of nuclear matter 
based on the nucleon interacting with the scalar a and the vector u meson fields linearly. 
And this is one of the models used in ref.0 to analyze (qq) p for p ~ 1.5 p a at T — 0. Note, 
however, that our results do not exactly match with theirs, since they use different values 
for the nuclear binding( -16 MeV and po — 0.17/m~ 3 ), but are in qualitative agreement. 

Spatially uniform nuclear matter, in spite of its simplicity, admits the following vital 
questions : 

(1) What happens to (qq) p when the a and uj fields are non-linear, inducing the 
couplings to be density dependent? These features were found to be very important say, 
in reducing the nuclear matter compressibility. 

(2) What is the effect on the effective nucleon mass? Is it similar to that of (qq) p ? 
What happens at T ^ ? 
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Keeping these in mind we choose, apart from Walecka, two other models used in H. 
So the models are : 

(1) Linear Walecka model where the coupling constants of the nucleon to o and u> 
fields, g a and g u respectively, remain constant with p. 

(2) ZM : (the usual one in the literature): It is constructed by changing the covariant 
derivative term in the Walecka model in such a way that, after an appropriate rescaling, 
the Lagrangian describes the motion of a baryon with an effective mass M* = m*M 
instead of the bare mass M. This information goes to the meson-baryon coupling, modify- 
ing it to an effective scalar coupling constant, dependent on p, while the vector coupling 
constant remains the same. 

(3) ZM3 : A variant of ZM is obtained if, instead of modifying the covariant derivative 
term, one simply modifies the kinetic energy term of the baryon. As before, after an 
appropriate rescaling, the Lagrangian describes a baryon of mass M*. This information 
manifests not only in the scalar-baryon coupling but also in the vector-baryon coupling. 
Both of them now depend on density. The vector and the scalar fields are now coupled. 

Connected in this way, 

£-ZM = £walecka(9a ~^ fl£) anC ^ £-ZM3 = C-W alecka{9 'a — > 9* a 9u ~^ fl£) 

The effective coupling constants are given by : g*/ 'g a — m* and g^/g^ = m* where 
M^/Mn = m*(cr) = (1 + g a o j M)~ x . The expression for the energy density at a given 
temperature T can be found in the mean field approach (MFA), 



where 



^ " (2k) 

Pb = j^t^ J d 3 k (n k - n k ). 



(2kY 

Here the degeneracy factor 7 = 4, n k and stand for the Fermi-Dirac distribution for 
baryons and antibaryons with arguments (E* =p v)/T respectively. The energy E*(k) = 
(k 2 + M* N 2 )2 and v is an effective chemical potential which preserves the number of baryons 
and antibaryons in the ensemble. The effective nucleon mass Mjy, is obtained from the 
minimization of £ . From eq. @ 

(qq) P ,T _ 1 m q d£ dM N | d£ dm^ | d£ dm a | d£ dg^ | d£ dg^_ 
(qq} m^U 2 dM N dm q dm^ dm q dm a dm q dg u dm q dg a dm q 

We adopt the rules [|], 

dm~ q = W N ^ and dm~ q = M^V (20) 
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and for the Walecka model the variation of the meson-couplings with m q is unspecified 
and neglected. Calculating the derivatives in eq . (|19"D , by using eq. (|T7D , we obtain a unified 
expression for (qq) at some T: 



= 1 ^[^W - M* N ) + (1 + a)^-(M N - M* N ) 2 - (1 + /?)-f^- 

(qq) m^U g* 2 g* 2 M N m 2 u M N 



(21) 

where, in terms of a and /3, the models are describes as : Walecka (a = (3 = 0), ZM (a = 
1 and/? = 0) and ZM3 (a = (3 = 1 ). The terms multiplied by a and f3 are obtained from 
the variation of the scalar and vector effective couplings with m q respectively. 



4. Results and Discussion for finite p 



The parameters of the models are presented in Table I to fix the binding energy per 
particle to be -15.75 MeV at p = 0.15/m -3 . It is found that the (qq) p /(qq) and cr e // 
are almost model independent upto a density p (Fig. 2). Throughout this paper a N 
= 45 MeV. With this, at p = p , a eff = 44.245 MeV for all the models. As expected 
0, the reduction is less than 2%. However, as evident in Fig.3, for high densities this 
reduction, 10 — 20% for ZM models, is not neglegible anymore. We stress here that 
(<?<?) Po l (qq)o = 0-69, for all the models. We calculate the contribution of 5£, using eq. (|i~6D , 
to the condensate and find it negligible at small densities. For example at p = p Q , 
(qq) S p £ I (qq) — 0.005 as compared to 0.69 for the total. Further, (qq) 6 p £ — — (AOMeV) 3 , 
reducing (qq) to — (203Mey) 3 from its vacuum value of -(230 MeV) 3 . At this density, 
the kinetic and the interaction energy almost cancel each other and there is saturation. 
This is reflected on (qq) also. The stabilizing density dependence of g* and make 

| (qq) | smaller than that in the Walecka model. Thus (qq) p /(qq) reduces with density 
to be at par with the leading order as shown in Fig. 2 for ZM models. This is more 
pronounced for ZM3 Model where (qq) increases slowly with p and the corresponding 
a e ff (Fig.3) is decreased. 

If there is scaling as in eq.(|6| or 0) or in any other form, the effective nucleon mass 
vanishes when the condensate vanishes. This is awkward for nuclear physicists since 
the meaning of nuclei or nuclear matter composed of zero mass particles is ill-defined. Is 
there some model where the condensate will decouple ? 

Indeed, as shown in Fig. 1-2, the ZM models decouple a decreasing condensate from 
the effective nucleon mass whereas in the Walecka model at p > po, (qq) p / (qq) goes up 
and M^/Mn — * 0. In-medium quark condensate is governed by cr e // ( eq.(|13"D) which in 
turn is the derivative of (S/p) and not of (eq.(0)). This explains why the (qq) and 
might have different density dependence and therefore the decoupling occurs. After 
all, at p = p , when £ / p — M N = — 15.75MeV , (qq) is found to be the same for all 
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hadronic models with very different M^(p ). 

To investigate the role of variation of meson-couplings with m q in the behaviour of 
(qq) p , we have separated the contributions of these variations in Fig. 4a and 4b respec- 
tively. It is interesting that in the ZM3 model the two contributions almost cancel each 
other upto p = 4p D which means that in eq. (^T|) the term actually important is the first 
one 

^(M N -M* N ) = ^a (22) 

which is proportional to the sigma field. In this type of model the a field is related not only 
with the scalar density (as in Walecka and ZM models) but have a new term depending 
on p, which comes from the new o~ — uj interaction present in this model J7| . So we can 
conclude that the condensate goes to zero in this model only because the contributions 
coming from the p\ are cancelled by the second term in eq. (^l|) which is proportional to 

2 ™ 2^2 



,*2 ; 



gfM iM N - M- nY = (23) 



For the ZM model the variation of the scalar coupling with m q is essential to get value of 
the condensate near the leading order behaviour. However, comparing the linear Walecka 
model with the non-linear models we use, we can conlcude that the effective density 
dependent coupling constants present in these models are mainly responsible for the con- 
densate not going up at intermediate densities (p/p = 1 — 2.5) as occurs in the Walecka 
model. 

In this context it is useful to mention that the density dependence of the coupling 
constants may be traced back to Brueckner Hartree Fock approaches [13| . In ZM models 
this dependence arises quite naturally from an interesting theoretical scaling property. It is 
noteworthy that in the original proposal of ZM3 [0], the authors claim that the justification 
for this model has to be understood in the spirit of chiral symmetry restoration. The 
present work confirms this, the (qq) being restored to zero at p ~ p Q . 

The non-linear interaction terms, in a heuristic way incorporates the effect of many- 
body forces leading to a suppression of the scalar field (Fig. 5). As a result, the lower bent 
of the nucleon mass, M/v, is checked in the medium. Moreover, in ZM3, even the vector 
field gets suppressed - holding back the increase of the condensate. Thus the condensate 
goes down at high density whereas the nucleon mass does not change so much. So, we 
have a working model of relativistic nuclear matter which has the right behaviour for (qq) 
built into it. Also it gives us the hope that one may indeed construct more refined models 
with the same features which combines the two desirable features of vanishing (qq) and 
non-vanishing effective nucleon mass. 
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5. Chiral Symmetry Restoration and Nuclear Matter Equation of State 

Let us analyse our results in more details. From Fig. 2 we find that for the ZM3 model 
which is nearest to leading order, the a e ff changes only 10% from at the critical density 
p c ~ 3.8po- This means the contribution from interaction 88 must remain small even for 
high densities. This can be seen clearly from the eq. ( |i~4"| ) when it is rewritten as 

djS/p) = _ d{SS/p) ^ _ (5S/p) 

dM N dM N M N 1 ' 

Indeed, this is the case for ZM models as shown by the equation of state (EOS, in 
short) given in [[/J, where S / ' p — M^ remains closer to zero even at p ~ 3 — 4p . In Walecka 
model, however, this is not true. S/p starts rising very sharply at p = 1.5po- Obviously, 
the vector repulsion takes over the scalar attraction in later. In ZM3, where both fields 
are supressed at high density, it is no wonder that the interaction part remains subdued. 

It is noteworthy that the relation given by the eq. (p4[) is exact for p = po- This 
explains the existence of the model independent point p where (qq) = 0.69 and a e ff = 
44.245 MeV, keeping the binding energy per particle fixed to -15.75 MeV and = 939 
MeV. 

^From the eq.(JZ4]) we understand that <7 e // and (qq) are directly connected to the 
EOS. Softer it is, the condensate moves closer to the leading order, realizing the restoration 
of chiral symmetry at high density. There is one parameter which governs the high density 
behaviour of the EOS for our models - that is the incompressibility, K = 9dP/dp where 
P is the pressure. Not only its value at po but also its variation with p should be small 
to obtain a softer equation of state. This is exactly what happens for ZM models, where 
incompressibility is small (K(po) = 156MeV for ZM3 and 220 MeV for ZM) and there 
is a possibility of chiral symmetry restoration. But this is impossible for Walecka model 
with its very stiff EOS. With these observations we can indeed connect the behaviour of 
the quark condensate with the nuclear matter parameters. 



6. Results and Discussion for finite T 



For finite temperature T, one must take into account the effect of pion gas which 
dominates the finite T vacuum fl4j| , |TJ]. This is easily done in the models considered 
above using the factorization hypothesis, assuming that the pion gas is non-interacting 
with the nucleons at high temperature. At moderate T this may introduce errors but 
at high T, near the chiral restoration point surely this is not a bad approximation. For 
moderate T careful calculations have been done by Leutwyler and Smilga Jl6| using the 
pion-nucleon scattering data through dispersion relations and recently by Koike |I7J using 
QCD sum rules including the effect of pion-nucleon scattering. These studies show that 
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the nucleoli mass does not change by more than a few percent at finite T, thus supporting 
our effective model calculations. 

In the present paper we restrict ourselves to the case of zero baryon density and finite 
T. This is the situation where the vector meson does not couple to the nucleon-antinucleon 
soup, but the scalar density is non-zero and an expression for the cr N (T) can be found 
easily. Thermal pions in the vacuum can be effectively taken through the T- dependence 

of/. H, 

fAT) = fAl - -^pf 1 ) (25) 

where 

B x {z) = - 2 dytf-z 2 )-*—-. (26) 
7r z Jz e y — 1 

The definition of the <Jn{T) for zero chemical potential is as follows : 

1-P,^L (27) 



<M>, ™rn„*f„ 2 (T) 



and 



p s a N (T) = i^(M N - M* (T)) + (1 + a)^-(M N - M* N {T)f}a N {T = 0). (28) 
where 

is the scalar density. 

Using M* N {T) = M N -g*a{T) and m* = M* N {T)/M N in eq.© where the scalar field 
cr(T) = p s m*" g*/m a 2 the sigma term reduces to a simple form : 

a N (T) = (2 + a - (1 + a) m*) m* a a N (0) , (30) 

In Figs.6(a-c), the change of quark condensate with T is shown, - with and without 
the pion effects in the vacuum. As in the case of finite density the behaviour of the 
quark condensate and the with temperature is found to be different. The condensate 
goes to its zero value in all the three models, - signalling perhaps, a new phase. The 
change in in ZM models is very little compared to that in Walecka model. In the 
latter, the condensate goes down very sharply, -just like the M^(T). The fall off is not 
compensated by the rise in p s : non-linearity in ZM models stabilizes <Jn{T) to 50 MeV 
at critical temperature (Fig. 7) which does not happen to linear Walecka model. This is 
also reflected in the a^{T) vs. T plot (Fig. 8). For this model ajq{T) increases sharply, 
meaning a sharp decrease in the condensate which indicates a first order phase transition 
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(Fig. 6a) at T c ~ 180 MeV. The ZM models, on the other hand, suggest a continuous 
and therefore a higher order transition. The inclusion of pions in the vacuum reduces the 
critical temperature (T c ) by 10%. For example in ZM3 model, T c becomes 230 MeV from 
its pionless vacuum-value T c = 260MeV. 

7. Summary and Conclusion 

In summary it is intriguing to find that the quark condensate at normal p is indepen- 
dent of the hadron-model chosen, - models with very different and incompressibility. 
But at high p and T it indeed depends on the modelling and it is possible to get a de- 
coupling of the condensate from the effective nucleon mass. The fact that {qq) p decreases 
with density in the ZM3 model (less than linearly) shows that EOS must be soft favour- 
ing a smaller incompressibilty (less than 220 MeV). However, the decrease in {qq) p does 
not compensate the linear rise of the density. Consequently, <r e // at high p reduces, by 
10 — 20%. The reduction at po is rather small in the hadron models we used. Further, pi- 
onic corrections arising from o and uo fields will be there- but as commented and observed 
in the reference - a calculation of a e ff with linear sigma model - their net effect may 
be small. 

With temperature, the condensate does go to zero. The ctjv(T) changes more dra- 
matically only for the Walecka model which shows a sharper increase compared to the 
non-linear hadron models where we have an increase of about 10% at T = T c , suggesting 
the latter might have higher order phase transition. 
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Figure Captions 

Fig. 1 (a-b) : Scaling laws of eq. (^j) and eq. (|7]) respectively for Walecka and Zimanyi- 
Moszkowski models (ZM and ZM3). 

Fig. 2 : Ratio of the condensates to the vacuum value as a function of density for 
three models are shown. 

Fig. 3 : a e ff with density. 

Fig. 4 : Leading order and other contributions to condensate ratios of fig. 2 for (a) : 
ZM and (b) : ZM3 models. 

Fig. 5 : The effective nucleon mass in the three nodels. 

Fig. 6 : A plot of the ratio of the condensate with the vacuum value and the effective 
mass with temperature for for (a) : Walecka, (b) ZM and (c) : ZM3 models. 
Fig. 7 : A plot of <jn with scalar density for three models. 
Fig. 8 : A plot of (7n with temperature for three models. 

Table 1. Coupling constants C%, C%, m*(p ) and (qq) / (qq) for different models 
with a fixed S/p - M N = -15.75 MeV at p a = 0.15/m" 3 . 



models 


cl 


cl 


m*(p ) 




Walecka 


352.69 


269.70 


0.54 


0.69 


ZM 


177.55 


63.47 


0.85 


0.69 


ZM3 


440.43 


303.22 


0.72 


0.69 
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